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ABSTRACT 

We present spectroscopic confirmation of two new lensed quasars via data obtained 
at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from 
the Dark Energy Survey (DES) and WISE based on their multi-band photometry and 
extended morphology in DES images. Images of DES J0115—5244 show two blue 
point sources at either side of a red galaxy. Our long-slit data confirm that both point 
sources are images of the same quasar at z s = 1.64. The Einstein Radius estimated 
from the DES images is 0.51". DES J2146—0047 is in the area of overlap between DES 
and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES 
and SDSS images. The SDSS fiber spectrum shows a quasar component at z s = 2.38 
and absorption compatible with Mg II and Fe II at Zi = 0.799, which we tentatively 
associate with the foreground lens galaxy. The long-slit Magellan spectra show that 
the blue components are resolved images of the same quasar. The Einstein Radius 
is 0.68" corresponding to an enclosed mass of 1.6 x 10 11 M 0 . Three other candidates 
were observed and rejected, two being low-redshift pairs of starburst galaxies, and 
one being a quasar behind a blue star. These first confirmation results provide an 
important empirical validation of the data-mining and model-based selection that is 
being applied to the entire DES dataset. 

Key words: gravitational lensing: strong - quasars: emission lines - methods: ob¬ 
servational - methods: statistical 
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1 INTRODUCTION 


Gravitationally lensed quasars provide unique insights into 
a variety of fundamental open problems in cosmology and 
extragalactic astrophysics (e.g. |Courb i n, Saha fe Schechter| 
20021. When a quasar is strongly lensed by a galaxy, it re¬ 


sults in multiple images of the same source, accompanied by 
arcs or rings that map the lensed host of the quasar. The 
light-curves of different images are offset by a measurable 


time-delay (e.g. Schechter et al. 1997 Tewes et al. 20131 
that depends on the cosmological distances to the lens and 


the source and the gravitational potential of the lens (Refs- 


dal 1964 1 . This enables one-step measurements of the ex¬ 


pansion history of the Universe and the dark matter halos 


of the massive galaxies that act as deflectors (e.g., Suyu 
et al. 120141. The microlensing effect on the multiple quasar 


images, induced by stars in the deflector, provides a quanti¬ 
tative handle on the stellar content of the lens galaxies (e.g., 
Schechter fe Wambsganss||2002| |Oguri, Rusu fe Falcol|2014| 


Schechter et al.|2014 Jimenez-Vicente et al. |2015[ > , and can 

simultaneously provide constraints on the inner structure 
of the lensed quasar, both the accretion disk size and the 


thermal profile (e.g. Poindexter, Morgan & Kochanek 2008 


Anguita et al. 2008 Eigenbrod et al. 2008 Motta et al. 2012 


as well as the geometry of the broad line region (e.g. Sluse 
et al. 2011, Guerras et al. 2013, Braibant et al. 2014). Fur¬ 
thermore, milli-lensing via the so-called flux ratio anomalies 
provides a unique probe of the mass function of substructure 


and thus ultimately of the nature of dark matter i 

Mao & 

Schneider 1998||Metcalf & Madam2001 Metcalf|2002 Dalai 

& Kochanek 

2002 

Nierenberg et al. 

2014). Finally 

source 


reconstruction of the lensed quasar and its host give a di¬ 
rect view of quasar-host coevolution up to s ~ 2 (Peng et al. 
2006 Rusu et al.|2015l. 


Advancement in the field is currently limited by the 
paucity of known systems suitable for detailed follow-up 
and analysis. A large sample of new systems will be trans¬ 
formative. To accomplish this, the STRIDES project]]] a 
broad external collaboration of the Dark Energy Surve)0 
(DES, http: //www. darkenergysurvey. org/index. shtml), 
aims at the discovery of the ~ 100 lensed quasars with pri¬ 
mary image brighter than i = 21 mag predicted by |Oguri fc| 
Marshall (20101 to lie within the DES footprint. 


The identification of such rare systems over 5000 square 
degrees of DES imaging data is a classic needle-in-a-haystack 
problem. The challenge is possibly greater than that faced 
in the SDSS dataset ( |Oguri et al.||2006| jlnada et al.||2008| ), 
as DES catalogs do not contain u-band data and there is 
no built-in spectroscopic dataset to aid in the quasar se¬ 
lection. Fortunately, the DES image quality (median seeing 
~ 0.9”) is better than that of SDSS, so one can rely on more 
accurate morphological information for candidate selection. 
New techniques have been developed in order to address this 
challenge ( Agnello et al.|2015 Chan et al.|2015 1. 

Here we report on the first spectroscopic confirmation 
of lensed quasars from DES. We have obtained spectra of 
five of the 68 high-grade, small-separation candidate lensed 


1 STRong lensing Insight in the Dark Energy Survey, PI Treu, full 
list of Co-PIs and Co-Is at http: //strides.physics .ucsb. edu 

2 See F laugher et al.| ( |2015| for the technical setup of the Dark 
Energy Camera. 


quasar s in the year-1 D ES data release footprint (hereafter 
Y1A1, Diehl et al.|2014 1, covering ~ 1200deg 2 in the South¬ 
ern Hemisphere. 

This Letter is organized as follows. Section [2] briefly il¬ 
lustrates the candidate-selection process and DES images of 
the two successful candidates from the first spectroscopic 
follow-up. Section [3] shows the long-slit spectroscopic data 
obtained for these systems. In Section[4]we present the lens¬ 
ing properties that can be inferred with current data and 
conclude in Section [5] Throughout this paper, DES and 
WISE magnitudes are in the AB and Vega system respec¬ 
tively. When needed in Sect. [4] we adopt a concordance cos¬ 
mology with f l m = 0.3, Ha = 0.7, Ho = 70kms -1 Mpc~ 1 . 


2 QUASAR LENS CANDIDATE SELECTION 


We selected small-separation candidates from the DES 
Y1A1 data release, using a combination of colour cuts and 
data mining and model-based selection. Following |Agnello| 


et al. (20151 we adopt a multistep strategy, whereby targets 


are first selected from catalog data, and candidates are then 
selected from the targets by analysis of the actual multi¬ 
band images. This multistep procedure allows one to keep 
the problem computationally fast. 


2.1 Preselection 


Objects are preselected in the DES catalogue based on their 
‘blue’ colours and extended morphology. The colour selec¬ 
tion, satisfying 


g — r < 0.6, r — i < 0.45, i — z < 0.55, 

2.5 < i — Wl < 5.5, 0.7 < Wl - W2 < 2.0, 

g — i < 1.2(* — Wl) — 2.8 , (1) 

is only used to exclude the majority of galaxies and nearby 
blue stars from the pool being examined. Of these, we retain 
those satisfying 

psf _mag — modeljnag ^ dmag (2) 


in DES g,r,i bands simultaneously, with dmag = 0.125 in 
g, r bands and 0.2 in i band. This ensures that the objects 
are not point-like (cf Reed et al. 2015, who use a similar but 
converse criterion for point-like objects). 

Out of ~ 2 x 10° ‘blue’ systems, extended in i —band and 
with acceptable colours in DES g, r, i, z and WISE Wl, W2 
magnitudes, ~ 4000 are brighter than limiting magnitudes 
i = 21 or Wl = 17 and extended also in g, r bands. 


2.2 Targets 

Of the 4000 blue extended objects, targets are selected by 
neural network classifiers (ANNs) based on their photome¬ 
try and multi-band morphology obtained at catalogue-level. 
In particular, DES model g, r, i, z magnitudes and WISE 
wlmpro (Wl), w2mpro (W 2) magnitudes are used for the 
photometry and DES model position angles and axis ratios 
in g, r, i, z are used for the morphology. ANNs produce mem¬ 
bership probabilities for each object to belong to different 
classes, from which the targets are selected with cuts in out¬ 
put probabilities. The full procedure is discussed in detail 
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Figure 1 . DES multi-band single-epoch images with best seeing of DES J0115—5244 (top) and DES J2146—0047 (bottom), in 5" X 5" 
cutouts. North is up and East is left. From left to right, DES g, r, i, z, Y. and colour-composites in gri, riz, izY. 
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Table 1 . Measured parameters of the systems in this paper. The redshift of source z s and lens z\ : are obtained from our long-slit spectra 
and the SDSS public data. The relative astrometry and magnitudes of brighter (A) and fainter (B) quasar image and lens galaxy (G) 
are obtained by modelling the single-epoch DES images with best image quality (quoted under ‘seeing FWHM’). The right ascension 
(resp. declination) increases to the E (resp. to the N). Upper limits are given for magnitudes with ^ 0.33 uncertainty. The galaxy in DES 
J2146—0047 is above the detection limit in DES Y band and in CFHT i. z bands. 


by Agnello et al. (20151. ANNs relying just on the photom¬ 
etry select 430 targets, out of which 136 are also selected by 
ANNs relying on multi-band morphologies as well. 


2.3 Candidates 

The DES g, r, i, z, Y image cutouts of selected targets have 
been modelled as combinations of point sources and galax¬ 
ies to extract the SEDs of the extended and point sources. 
Based on the model results, we ranked candidates from ‘0’ 
(not a lens) to ‘3’ (good lens candidate), based on whether: 
(i) the point sources have consistent SEDs at least in r, i, z 
bands; (ii) their SEDs are compatible with those of quasars; 
(iii) the fainter quasar image is also redder in g — r; (iv) 
a possible lens galaxy is detected in the residuals or in Y 
band. At this stage, the selection shifts from automatic pro¬ 
cedures (preselection, ANNs, modelling) to the investigator 
discretion, even though the ranking is based on the four 
quantitative criteria listed above. The final sample has 23 
grade ‘3’ candidates and 45, 51, 17 with grades ‘2’, T’, ‘0’ 
respectively. 


2.4 First two confirmed lenses 

Figure [l] shows DES multi-band, single-epoch images in 
the best seeing conditions of the two successful candi¬ 
dates for which we describe spectroscopic follow-up in the 
next Section. They are DES J0115—5244 at 01:15:57.32 
-52:44:23.20 (top) and DES J2146-0047 at 21:46:46.04 
-00:47:44.3 (bottom). 

The photometry and relative astrometry of the compo¬ 
nents in our composite models are summarized in Table [I] 
We have used the single-epoch DES cutouts with best image 
quality and adopted DES models of the PSF from nearby 
stars. For DES J2146—0047, we also fit Moffat profiles with 
the same structural parameters to each quasar image and a 
Sersic profile to the lensing galaxy in archival CFHT Mega- 
Cam i, z-band images (Programme 10BC22, PI: van Waer- 
beke). The resulting photometry is given in Table [l] below 
the DES one. Even though the galaxy is below the detec¬ 
tion limit in DES single-epoch data except in Y— band, it is 
securely detected in CFHT data in i, z bands as well. 

Three additional candidates have been observed during 
the same observing run and ruled out as lensed quasars. 
They are briefly described in Sect. |3.3] below. 
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Figure 2. ID sky-subtracted spectra of DES J0115—5244 (top) 
and DES J2146—0047 (bottom). The red and blue histograms 
show the spectra of the two components, demonstrating that they 
are images of the same quasar. Emission and absorption lines are 
marked by blue (solid) and red (dashed) vertical lines, respec¬ 
tively. 


3 SPECTROSCOPIC CONFIRMATION 


Spectra were obtained with the Inamori Magellan Areal 


Camera and Spectrograph (IMACS, Dressier et al. 20111 
on the Baade 6.5 m telescope on UT 2015 June 19. The f/4 
camera was used with 0.70” slit mask. The data were binned 
factors of 2 along the slit, giving a scale of 0.22”/pixel, and 
4 in the spectral direction, which combined with a 300 1/mm 
grating gave a dispersion of 2.92A /pixel. Two CCDs cov¬ 
ered the spectra over 3770 - 5260 A and 5350 - 6880 A. 
The seeing ranged between 0.8” and 1.2”, and all our candi¬ 
dates were partially or completely resolved. Two spectra of 
600s were taken for each object. The data were reduced us¬ 
ing standard IRAF routines. One-dimensional (ID) spectra 
have been extracted by modelling the 2D spectra as super¬ 
positions of Gaussian tracks in the spatial dimension, one 
per component, having peak positions (pi,P 2 ) that are lin¬ 
ear functions of the wavelength with the same slope (dp/dA). 
Even though the individual tracks are generally well sepa¬ 
rated, this procedure ensures that the resulting ID spectra 
are as independent as possible from one another and exploit 
all the information available in the 2D tracks. 

Five grade-‘3’ candidates were visible and observed. The 
resulting spectra are described below. The results on the two 
confirmed lenses are shown in Fig. [2] 


3.1 DES J0115—5244 

The ID spectra of DES ,10115—5244 show the same broad 
emission lines (C IV, He II, C III]) at z s = 1.64, with a uni¬ 
form ratio between the two spectra in the red, as shown in 
Figure [2] indicating that the low-order differences between 
the two spectra are due to differential reddening, or per¬ 
haps chromatic microlensing. Together with the presence of 
a red galaxy in the DES images, the spectra confirm DES 
J0115—5244 as a strongly lensed quasar, with image sepa¬ 
ration « 1.04”. Unfortunately, the S/N-ratio is too low to 
securely detect stellar absorption lines from the deflector. 


3.2 DES J2146—0047 

The ID spectra of DES J2146—0047 show two components 
with the same broad emission lines at z s = 2.38, consis¬ 
tent with the public SDSS fiber spectrum. Prominent Mg 
II and Fe II absorption lines at zi = 0.799 are detected in 
both spectra (fig. |2]). The ratio between the two spectra is 
constant, which together with the detection of a galaxy in 
izY bands confirms DES J2146—0047 as a strongly lensed 
quasar, with image separation ss 1.32”. We associate the 
lens redshift with the galaxy responsible for the absorption 
lines at zi. This system has been independently identified 
as part of the SDSS-III quasar lens sample (More et al., in 
prep). 


3.3 False Positives 


Two candidates, at 20:53:56.5 —56:09:36 and 22:17:52.5 
— 53:57:15 respectively, are pairs of compact, star form¬ 
ing galaxies. The third rejected candidate, at 22:00:24.11 
+01:10:37.56, is an alignment of a z = 1.37 red quasar and 
a blue star, with the same r — i and i — z colours by coinci¬ 
dence. 

In general, pairs of compact star forming galaxies at 
a ~0.2-0.3 are the main residual contaminant from the can¬ 
didate selection procedure, because their broad-band colours 
are consistent with those of quasars. This includes initial 
candidates that we rejected based on their available SDSS 
spectra. In WISE magnitudes, they tend to lie at W2 > 14 
and W1 — W2 < 0.8, across the limiting locus of |Assef| 
et al. (20131, a region that however is occupied also by 


spectroscopically-confirmed quasars and some of the can¬ 
didates from this search. 


4 LENSING PROPERTIES 


Even though the current data are not sufficient for a detailed 
lensing model, they can be used to derive simple proper¬ 
ties of the lens galaxies. To this aim, we adopt a Singular 
Isothermal Sphere (SIS) model for the mass density profile 
of the deflector, which is commonly considered the simplest 
model apt to describe galaxy-scale lenses ( Treu|[20 10). The 
SIS projected surface density is 


E (R) = -Y, cr (R/R E y 


(3) 


with Ecr = c 2 D s /(4nGDiDi s ) in terms of angular-diameter 
distances to the source ( D s ), to the lens (Di) and between 
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DES J0115—5244 

DES J2146—0047 

Re 

0.513" 

0.680" 

M tot,i 

26.0T20.0 

8.4T3.2 

A^totjC 

16.0T13.0 

6.6T1.4 

tTais (km/s) 

190 

215 

Me (M 0 ) 

0.95 x 10 11 

1.6 x 10 11 


Table 2. Estimated lensing parameters and properties of the lens 
galaxy: Einstein radius Re (in angular units) and total magnifi¬ 
cation fitot,i from the fit to the quasar image positions and fluxes; 
magnification /xtot,c on the quasar continua, from the ID spectra; 
and SIS velocity dipersion cr s i s and enclosed mass Mg from the 
measured redshifts. For DES J0115—5244, cr s i s and Me are un¬ 
der the hypothesis that the lens is midway between observer and 
source (i.e. zi = 0.635). 


lens and source ( Di s ), and the Einstein Radius Re is defined 
such that it encloses a projected mass 

Me = 7rE cr R| . (4) 


The expected line-of-sight velocity dispersion of stars in the 
lens galaxy can be estimated via 


_ c 2 Re D s 
An DiD i 


Finally, neglecting microlensing, the total magnification of 
the point source can be estimated from the flux ratio fs / /a 
as 

_ n 1 + Sb/U ( 6 \ 

Mtot —r< 2 1 — f g /f A ’ (6) 

where A is the brighter quasar image. Then, from Table [l] 
and confirmation spectra, we can estimate Re, Me, <r s is and 
AHot- The Einstein radius is approximated by half the image 
separation (an exact result for SIS profiles). For the flux ra¬ 
tio, we use the average between i and 2 band, for which the 
photometric fit is more robust against PSF mismatch and 
signal-to-noise ratio, and the effects of dust and microlens¬ 
ing are minimized. The results are summarized in Table [2] 
For DES J0115—5244, we adopted a lens redshift zi = 0.635, 
which minimizes the critical surface density E cr and yields a 
lower bound on <r a i s . We also list the magnification inferred 
from the ratio of continua in the spectra around 6000 A. 
With the positional uncertainties from Table [T] the relative 
uncertainty on Re is 1%. With current data, the uncertain¬ 
ties on the magnifications are of the same order as their 
average value. 


5 CONCLUSIONS 

We have presented the very first results of a campaign to ex¬ 
tend the known samples of lensed quasars by exploiting the 
large footprint of wide-field photometric surveys. In partic¬ 
ular, we have spectroscopically confirmed two new systems, 
found in the DES Y1A1 release with data mining techniques. 
Two additional candidates were ruled out by spectroscopic 
observations as pairs of z « 0.2 compact, narrow-line galax¬ 
ies. Another candidate was ruled out as an alignment of a 
quasar and a blue star. The basic characteristics of the two 
confirmed systems are as follows. 

DES J0115—5244 Consists of two images of a quasar 


at z s = 1.64 lensed by a foreground galaxy visible in the 
DES images. No redshift is available for the lens galaxy. 
The Einstein Radius is estimated to be 0.51”. 

DES J2146—0047 Consists of two images of a quasar 
at z s = 2.38, with prominent Mg II absorption at zi = 0.799, 
which we tentatively associate with the lens galaxy redshift. 
The redshift of the lensing galaxy is not confirmed in the 
spectroscopic follow-up by More et al. (in prep), since detec¬ 
tion of the Mg II and Fe II lines depends on S/N ratio and in 
fact is less evident in the worse of our two exposures, where 
track deconvolution becomes more noisy. Thus, deeper spec¬ 
troscopic data will be needed to measure the redshift of the 
deflector, possibly based on stellar absorption lines. The Ein¬ 
stein Radius is estimated to be 0.68”. 

The main class of contaminants (including candidates 
rejected based on their SDSS spectra) consists of groups of 
compact, star-forming galaxies at z ~ 0.2 — 0.3, because 
of their broad-band colours and compact morphology. With 
the observed sample, a strict cut in WISE W1 — W2 vs W2 
would give a 100% success rate, which drops at 40% when 
the cut is relaxed. Still, with just five systems currently at 
hand, we would rather caution against these simple esti¬ 
mates. 

Our search has delivered over 100 additional candidates 
from the Y1A1 data - and we expect many more from the 
next seasons of DES data. The results of this first follow¬ 
up effort are encouraging. However, a systematic follow-up 
campaign is needed to confirm large numbers of candidates, 
assess the purity of our selection technique, and carry out 
the many scientific investigations enabled by lensed quasars. 
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